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! Abstract 

00 

I We have calculated the strength distributions of the giant monopole resonance in the even- A tin 

' isotopes {A = 112 — 124) which were recently measured in inelastic a-scattering. The calculations 

OO ■ 

, were performed within two microscopic models: the quasiparticle random phase approximation 

' (QRPA) and the quasiparticle time blocking approximation which is an extension of the QRPA 

including quasiparticle-phonon coupling. We used a self-consistent calculational scheme based on 



. the HF+BCS approximation. The single-particle continuum was exactly included on the RPA 

level. The self-consistent mean field and the effective interaction were derived from the Skyrme 
energy functional. In the calculations, two Skyrme force parametrizations were used. The T5 
parametrization with comparatively low value of the incompressibility of infinite nuclear matter 
(i^oo = 202 MeV) gives theoretical results in good agreement with the experimental data including 
the resonance widths. 

PACS numbers: 21.60.-n, 24.30. Cz, 25.55.Ci, 27.60. +j 
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I. INTRODUCTION 



The investigation of the isoscalar giant monopole resonance (ISGMR), the so-called 
breathing mode, is one of the fundamental problems of nuclear physics. The energy of 
the ISGMR enables one to determine parameters characterizing the incompressibility of in- 
finite nuclear matter (INM), in particular, the value of K^o- These collective resonances 
can be studied experimentally in inelastic a-scattering at small angles (see, e.g., Ref. |l| 
and references therein). Theoretical investigations of these states are based mainly on the 
self-consistent microscopic approaches (see, e.g., Refs. nj, \m) 

including, first of all, scaling and constrained Hartree-Fock (HF) models and the random 
phase approximation (RPA) and on the Landau-Migdal approach that starts with a phe- 
nomenological single-particle basis and with the independently parametrized particle-hole 
zero-range interaction (see, e.g., Refs. Q,Q,Q and references therein). It is important 
to note that the incompressibility can not be measured directly but it can be deduced 
theoretically by comparing the experimental energies of the ISGMR with the correspond- 
ing calculated values. The most widely used approach is based on the self-consistent HF 
or RPA calculations of the mean energies of the ISGMR using effective Skyrme or Gogny 
forces. Because Koo can be calculated from the known parameters of the given force, its 
value is estimated as the one corresponding to the force that gives the best description of 
the experimental data. The non-relativistic estimates obtained in such a way lead to the 
value i^oo = 210 ± 30 MeV (see, e.g., Refs. |iQ,fly,0,|8|,|9|,Q), though the recent results 



TT . 

favor the upper limit of this estimate (see [ll|, ll2|). In the Landau-Migdal approach one 
obtains K^o from the scalar-isoscalar Landau-Migdal parameter /q. Here K^o was always of 
the order of 240 MeV Q- 

Note that within the relativistic mean-field (RMF) theory the INM incompressibility is 
usually restricted to the interval K^o = 260 ±10 MeV (see, e.g., Ref. 13]) that is considerably 
higher than the non-relativistic limits. However, recently it was obtained in Ref. that a 
zero-range (point-coupling) representation of the effective nuclear interactions in the RMF 
framework leads to the reduction of up to the value of 230 MeV. 

n 

In the present paper we investigate theoretically the new experimental data [IS] on the 
strength distributions of the ISGMR in the even-A tin isotopes (A = 112 — 124). This is the 
main goal of our work. The calculations are performed within the framework of the recently 
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developed microscopic model that takes into account the effects of the quasiparticle-phonon 
coupling (QPC) in addition to the usual correlations included in the conventional RPA. 

The paper is organized as follows. In Sec. [TTl the model is described. The particular 
attention is paid to the dynamical pairing effects which are important to solve the problem 
of the 0+ spurious state in the ISGMR calculations in open-shell nuclei. In Sec. IIIII we 
describe the details of our calculational scheme and present the results and their discussion. 
Conclusions are drawn in the last section. Appendix |X] contains auxiliary formulas. 



II. THE MODEL 



A. General scheme 



Two microscopic models were used in our calculations. The first one is the well-known 
quasiparticle RPA (QRPA). The basic ingredients of this approximation are the nuclear mean 
field (including the pairing field operator) and the residual particle-hole (ph) interaction. In 
the self-consistent QRPA these ingredients are related to each other by the consistency 
condition. The nuclei excitations are treated as superpositions of the two-quasiparticle (2q) 
configurations. This model is applicable to a wide range of nuclei including open-shell 
ones as the pairing correlations of nucleons are taken into account. The QRPA reproduces 
well the centroid energies and total strengths of giant multipole resonances but not their 
widths. In order to reproduce the total widths of the resonances it is necessary to enlarge 
the configuration space by adding 4q configurations, i.e. to extend the (Q)RPA. The most 
successful approaches in this direction are the models which take into account the QPC in 
addition to the correlations included in the (Q)RPA (see Refs. [l5|, [l9|, 120] and references 
therein) . 

In the present investigation the QPC contributions are included within the framework 
of the recently developed quasiparticle time blocking approximation (QTBA) which is an 
extension of the QRPA in this sense. On the other hand, since in the QTBA the pairing 



correlations are also inc^ 
decoupling of diagrams 

n 



uded, this model is a generalization of the method of chronological 



211 ] which is a base of the Extended 



'heory of Finite Fermi Systems 



23|. The basic equation of our 



[iq . Details of the QTBA model are described in Refs. 22, 
approach (both in the QRPA and in the QTBA) is the equation for the effective response 
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function R (ui). In the shorthand notations it reads (we will follow notations of Ref. [23|) 



R'^iu) = A{uj)- A{iu) T R'^'icu) (1) 

where A^u) is a correlated propagator and is an amplitude of the effective residual in- 
teraction. In the case of the QRPA, A{uj) reduces to the uncorrelated 2q propagator A{uj). 
In the general case including pairing correlations, the amplitude can be represented as a 
sum of two terms 

jr = jr(ph) + jr(pp) (2) 

where the amplitude J^^p*^) represents interaction in the ph channel and JF^pp) includes con- 
tributions of the interaction both in the particle-particle (pp) and in the hole-hole (hh) 
channels (in the following for brevity we will use the unified term pp channel implying also 
the hh-channel contributions). 



Let us emphasize that the general formulas of the QTBA derived in Ref. [22] are valid 
both in the self-consistent and in the non-self-consistent approaches. In the present paper, 
we use self-consistent calculational scheme based on the HF and Bardeen-Cooper-Schrieffer 
(BCS) approximations (in what follows, we will refer to this scheme as the HF+BCS approx- 
imation). The self-consistent mean field and the effective residual interaction were derived 
from the Skyrme energy functional by means of the known variational equations. In the 



calculations, the T5 and T6 Skyrme forces (see Ref. 2^) were used 



An important property of these parametrizations is that they produce the nucleon effec- 
tive mass m* equal to the bare nucleon mass m. This is a consequence of the fact that the 
T5 and T6 Skyrme- force parameters are constrained by the relations (see, e.g., 25^]): 

t2 = -7Tii 5 + 4x1 , X2 = --—-^. 3 

6 b + Axi 

In this case the contribution of the velocity-dependent terms (except for the spin-orbital 
ones) into the energy functional and the mean field reduces to the derivatives of the nucleon 
density, i.e. to the simple surface terms. As a consequence, contribution of these terms into 
the effective interaction derived from such an energy functional also has very simple form. 
To see this, consider the energy density Ti of the Skyrme energy functional S defined as 

S = [ drn{r). (4) 



4 



In the sufficiently general case it is given, e.g., in Ref. 2J]. However, if the equations ([3]) 
hold, the energy density acquires the form 

+ ±t, {3(1 + 1 xi)(Vp)2 - 3(xi + i) [(Vp„)2 + iVp.f] - xi J2 + Jl + Jl} 

- + \x2){Vpf + (X2 + \)[{Vpnf + {Vp,f] +X2J' + JI + jj} 

+ ii'^3[(l + |a;3)p'-(x3 + |)(p^ + pj)]p" 

+ i w^o ( J ■ vp + j„ ■ vp„ + jp ■ vpp) + T^^^^j + n^^.^ (5) 

where '^^q^^^j is the Coulomb energy density including the exchange part in the Slater ap- 



«c„„,W-^/^-'^^#^-?f!V^VfW, (6) 



proximation, i.e. 

Coui ^ ' 2, J |7» /''I 4\^7r 

and '^^pg^jj. is the density of the pairing energy. In the applications of the models based on 
the Skyrme energy functionals it is frequently taken in the simplest form 

^pair = i ^0 (x:X„ + (7) 

which was also used in our calculations. In Eqs. (l5])-(l7l), p = Pn + Pp, Pq, Tq, and Jg are the 
normal densities and Xg is the anomalous local density of the nucleons of the type q = n,p 
(neutrons or protons). In particular, pg is the local particle density, Tq is the kinetic-e nerg y 



density, and Jq is the spin density. They are defined in the usual way (see, e.g., Ref. |26| ) . 
In the case of the spherically symmetric nucleus and within the HF+BCS approximation 
they have the form 

p.(^) = (8) 

(1) 



(1) 



(1) 



2j _|_ X 

Yl K^i^i:^- "(1) ^(1) ^(1) • (11) 



(1) 



Here and in the following we use the notations of Refs. 22|, |23| for the single-quasiparticle 



basis functions in the doubled space ip-^ which are labelled by the composite indices 



1 = {[l],m]^} where [1] = {{l),!]^}, (1) = {qi,ni,li, j^}, and r]^ = ±1 is the sign of the 
quasiparticle energy E-^ = rj^Ef^^y I.e., the symbol (1) stands for the set of the single-particle 
quantum numbers excepting the projection of the total angular momentum m^, R^^^{r) is 



the radial part of the single-particle wave function, is the occupation probability, and 



As can be seen from Eqs. (jlj) and ([5]), the following equality holds 

6S !f 



constant . (12) 



In particular, it means that the equations of motion derived from such an energy functional 
S contain the nucleon effective mass m*{r) = m. The spin-scalar part of the effective 
interaction in the ph channel corresponding to S is determined by the relation 

This ansatz completely includes velocity-dependent contributions because of Eq. (fT2l) . In 
the explicit form we have 



(^i to (1 - xo) + ^ t3 |(1 + I xs) (1 + a) (2 + a) 

- (X3 + i) [2p" + 4ap„p"-i + a{a - l){pl + pJ)p"-2] |^ 6{r - r') 

+ ^ [h (1 + X2) - ti (1 - xi)] A 5{r - r') , (14) 

^o^piry) = (^,{l + \x,) + l^h\^{l + \x,){l + a){2 + a)p- 

- (X3 + \)a [{a + l)p" - 2{a - l)p„ p^ p""^] |^ 5{r - r') 

+ I [t2 (1 + \x2) - 3ti (1 + \xi)\ A 5{r - r') . (15) 

The formulas for the components ^^^l and J^l^^^ obtained from Eqs. ffT^ and f[T^ by 
replacing indices n by p and p by n and by adding the Coulomb interaction to J-'o^pp- 

Let us note that in addition to the simplicity of the formulas for the effective residual 
interaction there exist the physical reasons to use the Skyrme forces with m* = m. It 
is known that for heavy and medium mass nuclei the single-particle spectra obtained in 
the HP calculations with such Skyrme forces better reproduce the experimental energies as 
compared with the case of the forces with m* /m ~ 0.7. This results in better description 



of the excitations of the even-mass nuclei in the RPA and QRPA. The same is true for the 
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23|) is used. 



QTBA if the subtraction procedure (see Eq. (|T9|) below and Refs. 

The spin-vector components of the effective interaction are not determined uniquely from 
Eq. which is valid only for the spin-saturated nuclei. In our calculations these components 
are taken in the simple form of the Landau-Migdal zero-range force with known parameters 
Co, g, and g' (see, e.g., Ref. 3]). However, the spin-vector components of the interaction 
do not enter equations for the O"*" excitations. 



Note that in the fully self-consistent HF-RPA calculations 27|] it was obtained that the 
spin-orbital and the Coulomb components of the effective residual interaction slightly shift 
the mean energy of the ISGMR in the opposite directions. This effect is also confirmed 
by our estimates of the mean energies of the ISGMR based on the constrained HP method 

n n 

described in Ref. |3[ (see also [10| for details of our calculational scheme). Namely, in the 
test calculations for ^^°Sn we have obtained the following results using T5 Skyrme force and 
neglecting pairing correlations. The constrained HP method (fully including the spin-orbital 
and the Coulomb contributions) yields for the mean energy A/mi/m_i the value 15.1 MeV 
where mi and m_i are the energy-weighted moments defined by Eq. ( 122|) below for the 
infinite energy interval {Ei = 0, E2 = 00). For the same interval we obtain in the RPA the 
value 15.4 MeV if we neglect only the spin-orbital contributions in the residual interaction. 
If we neglect both the spin-orbital and the Coulomb contributions in the RPA interaction 
we obtain the value 15.0 MeV that differs from the constrained HP result only by 0.1 MeV. 

In our QRPA and QTBA calculations, the spin-orbital components of the effective inter- 
action (but not of the mean field) are neglected. Por the reasons described above we also 
neglect the Coulomb contribution into J^lfpp- The effective interaction in the pp channel 
and the gap equation within the HP+BCS approximation are determined by the formulas 
of Appendix A of Ref. [23| with jF5(r) = |Vo (see also Appendix |A] of the present paper). 



B. Dynamical pairing effects in QRPA and QTBA 

One of the important questions arising in the QRPA and QTBA calculations is the 
question of completeness of the configuration space. The size of the basis in this space has 
an impact practically on all the calculated quantities. In particular, configurations with 
a particle in the continuum are responsible for the formation of the escape widths of the 
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resonances. The well-known method to include these configurations on the RPA leve 



is the 



use of the coordinate representation within the Green function formalism (see Ref. [28|). We 
use this method in our approach as described in Ref. j^]- However, incorporation of the 
pp-channel contributions in the coordinate representation leads to considerable numerical 
difficulties. At the same time, the pp-channel contributions (so-called dynamical pairing 
effects) are very important in the calculations of O"*" excitations in the open-shell nuclei, first 
of all because of the problem of the O"*" spurious state. For this reason we have developed 
a combined method which is a modification of the so-called (r. A) representation proposed 



in Ref. [29|] for the QRPA problem. Within this method only the ph channel is treated in 
the coordinate space, while the dynamical pairing effects are included in the discrete basis 
representation. 

Considering the general case of the QTBA, note that taking into account decomposition 
([2]) one can rewrite Eq. ([1]) in the form 

R^^{uj) = A^'^^'^+PP) (u) - A^'^'+PP'^ (uj) J^(P^) R'^iuj) (16) 

where propagator 

^(res+pp) 

(cu) is a solution of the equation 

^(res+pp) ^ j^^^^ _ ^^^^ jripp) ^(res+pp) _ (^^j^ 

In the present work we use the version of the QTBA in which the ground state correlations 
caused by the QPC are neglected. In this case the correlated propagator A{u) is defined by 
the equation 

A{uj) = A{uj) - A{u}) $(cu) A{uj) (18) 
where A{uj) is the uncorrelated QRPA propagator, 

$(cj) = $('^<=^)(cj) -$('^<=^)(0), (19) 

and $*^'''^'^^(u;) is a resonant part of the interaction amplitude responsible for the QPC in our 



model (see Refs. 



22 



231] for details). Combination of Eqs. (1171) and (1181) leads to the new 



equation for A'^'^'^^^''^p\uj): 

^(rcs+pp) _ ^^^^ _ ^^^^ ^^^^ ^ _^(pp) j ^(res+pp) (^^^ _ (2O) 

As a result we obtain that the pp-channel contributions can be included by modification 
of the equation for the correlated propagator, i.e. by replacing Eq. (fT8|) by Eq. (!20|) . The 
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modification is reduced to the additional term JF(pp) added to the amphtude ^{u). The 
respective equations in terms of the reduced matrix elements are drawn in Appendix El It 
is worth noting that the QPC in the QTBA is included both in the ph channel and in the 
pp channel because there is no difference between these channels in the representation of 
the single-quasiparticle basis functions in the doubled space (t/'i, see Ref. ^]) which is used 
in Eqs. ( flTl) . (fTSi) . and ( !20l) . This is true both for the system of equations ([T]), (fTSl) and for 
the system (fTB]), (|2U1). 

Note, however, that in practice Eq. ( fT6|) for W^^lu) is solved in the coordinate represen- 
tation (to take into account the single-particle continuum), while Eq. fl20|) is solved in the 
restricted discrete basis representation. This fact greatly simplifies the problem as compared 
with the initial Eq. ([1]) in which both the ph-channel contribution and the pp-channel one 
are included in the coordinate representation. On the other hand, the use of the restricted 
discrete basis representation for the pp channel is fully consistent with BCS approximation 
in which the gap equation is solved in the same restricted basis. 

The general scheme described above ensures that the energy of the 0"*" spurious state 
(so-called ghost state) is equal to zero both in the QRPA and in the QTBA. However, there 
still remains the following problem: in the QTBA the ghost state can be fragmented due 
to its coupling to the 2q(S)phonon configurations, despite the energy of the dominant ghost 
state is equal to zero. It can lead to the spurious states at low energies distorting respective 
strength functions. In particular, these fragmented spurious states will produce non-zero 
response to the particle-number operator which has to be exactly equal to zero in a correct 
theory (as, for instance, in the QRPA including pp channel that was proved by Migdal, see 
jsojl). In the present calculations this problem is solved with the help of special projection 
technique which will be described in a forthcoming publication. 

III. CALCULATIONS OF THE GIANT MONOPOLE RESONANCE IN THE TIN 
ISOTOPES 

A. Numerical details 

The method described above has been applied to calculate the strength distributions of 
the isoscalar giant monopole resonance in the even-A tin isotopes {A = 112 — 124) which were 
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recently measured experimentally with inelastic scattering of a particles (see Ref. ISj). The 

3F+BCS approximation using 
24| including the pairing-force 



ground state properties of these nuclei were calculated within 
T5 and T6 Skyrme forces with the parameters drawn in Ref. 
strength Vq = —210 MeV-fm^ in Eq. ([7]). For all tin isotopes under consideration, the pairing 
window for the neutrons contains 22 states including all the discrete states and one or two 
quasidiscrete states. The criterion to select quasidiscrete states is described in Ref. |23|. 

To calculate the strength function of the ISGMR, the equation ([1]) for the effective re- 
sponse function B^'^iuj) was solved. The strength function S{E) is determined by Bf^iuS) 
via the formula 

S{E) = i-Im5^(er°)*,<3,(E + .A)(eV°),3 (21) 

1234 

where E is an excitation energy, A is a smearing parameter, is an external field, and 
e is an effective charge operator. In the case of the isoscalar 0^ excitations the one-body 
operator eV^ is proportional to the identity matrices both in the spin and in the isospin 
indices. Its radial dependence is taken in our calculations in the form eV^ = r^. The 
smearing parameter was taken to be equal to 500 keV that approximately corresponds to 
the experimental resolution for the data presented in Ref. ISj. 

In the calculation of the QTBA correlated propagator A{uj) entering Eq. ([T]), the valence 
zone for the neutrons coincides with the pairing window. The valence zone for the protons 
contains 20 states including all the discrete states and several quasidiscrete states as de- 
scribed in Ref. [23j]. Let us emphasize that the restricted valence zone is used only in the 
calculation of the discrete part of the propagator A{uj) including QPC effects and in the 
calculation of the phonons (see below). In the ISGMR calculations, the configurations with 
the particle in the continuum are included completely in the RPA-like part of A{uj) (see 
Ref. y for details). 

The set of phonons in the QTBA calculations included collective modes with values of 
the spin L in the interval 2 ^ L ^ 9 and with natural parity tt = (—1)^. The phonon 
characteristics were calculated within the QRPA using configuration space restricted by the 
valence zone described above. The maximal energy of the phonon was adopted to be equal 
to the value 10 MeV which is approximately equal to the nucleon separation energy for the 
given tin isotopes. The second criterion to include the phonon into the phonon space was 
its reduced transition probability B[EL) which should be more than 10% of the maximal 
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B{EL) for the given spin. According to these criterions, the total number of phonons 
included in the QTBA calculations is equal to 21 for ^^^Sn, 19 for ^^^Sn, 23 for ^^^Sn, 26 for 
ii^Sn, 29 for ^^osn, 27 for ^^^Sn, and 31 for ^^^Sn. 

To describe correctly effects of a fragmentation of the resonances in the QTBA arising due 
to the QPC it is very important to use the phonon space with the phonon characteristics 
close to the experimental ones. However, both the T5 and the T6 Skyrme forces do not 
provide satisfactory description of the experimental energies and transition probabilities 
within the self-consistent QRPA scheme presented in Section [III For this reason, in the 
calculation of the phonons (and only in this calculation) we have used the QRPA scheme 
which is self-consistent only on the mean-field level. More specifically, the mean field was 
calculated within HF-I-BCS approximation based on the T5 Skyrme force, while the effective 
residual interaction was taken in the form of the Landau-Migdal zero-range force with the 



standard set of the parameters (see, e.g., Ref. [31|), except for the parameter f^^. This 
parameter was adjusted for the each nucleus to reproduce the experimental energies of the 
2f and 3^ levels. As a result, the parameter f^^ takes the values in the interval —1.54 ±0.11 
for the phonons with the positive parity and the values in the interval —1.83 ± 0.06 for the 
phonons with the negative parity. 



B. Results and discussion 

The results for the ISGMR strength distributions in the even-A ^^^~^^^Sn isotopes are 
presented in Fig. [T]and in Tabled The energy- weighted moments were determined as 



rrik 



/ E^S{E)dE. (22) 

J El 



The energy interval limited by Ei = 10.5 MeV and E2 = 20.5 MeV was taken the same as 
in Ref. (l8|. As can be seen from the Tabled the agreement of the theoretical results with 
the experimental mean energies is good both in the QRPA and in the QTBA. The fact that 
the mean energies obtained in the QRPA and in the QTBA are very close to each other is 



explained by the subtraction procedure used in our calculations (see Eq. ( |T9|) and Ref. [23 1 
for the discussion). The main reason of the agreement with the experiment is that in this 
calculation we used the self-consistent scheme based on the T5 Skyrme-force parametrization 
with comparatively low value of the incompressibility of INM {Kao = 202 MeV). The other 
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10 20 
Energy (MeV) 



FIG. 1: Isoscalar giant monopole resonance in the even-A ^^^~^^^Sn isotopes calculated within 
QRPA (dashed line) and QTBA (solid line). The results are obtained within self-consistent 
HF+BCS approach based on the T5 Skyrme force. The smearing parameter A is equal to 500 
keV. Experimental data (solid squares) are taken from Ref. 181 ]. 
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parametrizations with K^o around 240 MeV give too large mean energies of the ISGMR in 
the considered tin isotopes as compared with the experiment. 

For comparison, in Table [III we draw the QRPA results obtained with the T6 Skyrme 
force {Koo = 236 MeV). As can be seen, the T6 mean energies mi /mo are greater than the 
experimental values for the tin isotopes by more than one MeV. This fact agrees with the 
results of Ref. 32j where it was obtained that the relativistic RPA calculations based on the 



force with K^o = 230 MeV consistently overestimate the centroid energies of the ISGMR in 
the same tin isotopes. 

Note that the value K^o = 202 MeV corresponding to the T5 Skyrme force lies within the 
interval 210 ± 30 MeV which was considered for a long time as the non-relativistic estimate 



for this quantity. The recent results 
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12l | giving Koo around 230-240 MeV were obtained 
on the base of the experimental data in fact only for the one nucleus ^°^Pb which is doubly 
magic. However, the question which arises is whether the doubly magic nucleus is the best 
candidate to determine the value of the INM incompressibility in view of the strong shell 
effects taking place in this case. On the other hand, by comparing the RPA results for ^ospj-^ 
shown in the Tables [T] and [IT] one can see that the T6 Skyrme force with Koo = 236 MeV 



nice 

see 



y reproduces the experimental data for this nucleus ((?^i/?^o)exp = 14.2 ± 0.3 MeV, 
l|), while the T5 force gives the result which is lesser by 1.3 MeV. Thus, the question 
about the precise value of Koo is not resolved within the framework of our approach. 

The theoretical values of the ISGMR widths F were obtained from the Lorentzian fit of 
the calculated functions S{E). In contrast to the mean energies, the QRPA and the QTBA 
give substantially different results for the width. It is well known that the spreading width 
F^ is a considerable part of the total width of the giant resonance. The QRPA does not 
produce F^-, while in the QTBA it is formed by the 2q(8)phonon configurations. This is the 
reason why the QRPA strongly underestimates the experimental values of F, while very 
good agreement is achieved in the QTBA. 

To investigate the nature of dependence of the ISGMR mean energies on the neutron 
excess {N — Z) we calculated the unperturbed 0^ response substituting the (Q)RPA uncor- 
related propagator A{uj) in Eq. fl2T|) instead of R^^{uj). This response corresponds to the 
independent quasiparticle model (IQM). The results are presented in Table UTT] in compari- 
son with the (Q)RPA results obtained in the same energy interval 10-30 MeV. This interval 
was chosen to exclude contribution of the low-lying strength arising in the IQM response. 
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TABLE I: Mean energies and widths for the ISGMR strength distributions in the even-^ ^^^~^^^Sn 
isotopes calculated for 10.5-20.5 MeV energy interval. Theoretical results are obtained within 
self-consistent HF+BCS approach based on the T5 Skyrme force (-fCoo = 202 MeV). Experimental 



181 ]. The RPA results for ''"'^ Sn and ^OSp^ ^j-g drawn for comparison. 



values are taken from Ref. 
The values for 208p|-, ^^^.g calculated for 5-25 MeV energy interval. 







\/ m-i /m„i , MeV 


mi /mo, MeV 


\/^m^x|m^ , MeV 


r, MeV 


lOOgj^ 


RPA 


16.4 


16.5 


16.7 


2.2 




OR PA 


16.0 


16.1 


16.3 


1.9 


"2Sn 


QTBA 


15.9 


16.0 


16.4 


3.9 




Experiment 


16.1 ± 0.1 


16.2 ± 0.1 


16.7 ± 0.2 


4.0 ± 0.4 




OR PA 


15.9 


15.9 


16.2 


2.0 


ii^Sn 


QTBA 


15.7 


15.9 


16.3 


3.9 




Experiment 


15.9 ± 0.1 


16.1 ± 0.1 


16.5 ± 0.2 


4.1 ± 0.4 




OR PA 


15.7 


15.8 


16.1 


2.1 


ii^Sn 


QTBA 


15.6 


15.7 


16.1 


4.0 




Experiment 


15.7 ± 0.1 


15.8 ± 0.1 


16.3 ± 0.2 


4.1 ± 0.3 




OR PA 


15.6 


15.7 


16.0 


2.3 


"8Sn 


QTBA 


15.5 


15.6 


16.0 


4.2 




Experiment 


15.6 ± 0.1 


15.8 ± 0.1 


16.3 ± 0.1 


4.3 ± 0.4 




OR PA 


15.4 


15.5 


15.8 


2.4 


i20Sn 


QTBA 


15.3 


15.5 


15.9 


4.3 




Experiment 


15.5 ± 0.1 


15.7 ± 0.1 


16.2 ± 0.2 


4.9 ± 0.5 




QRPA 


15.3 


15.4 


15.7 


2.5 


i22Sn 


QTBA 


15.2 


15.4 


15.8 


4.4 




Experiment 


15.2 ± 0.1 


15.4 ± 0.1 


15.9 ± 0.2 


4.4 ± 0.4 




QRPA 


15.2 


15.3 


15.6 


2.6 


i24Sn 


QTBA 


15.1 


15.3 


15.7 


4.4 




Experiment 


15.1 ± 0.1 


15.3 ± 0.1 


15.8 ± 0.1 


4.5 ± 0.5 


i32Sn 


RPA 


14.8 


14.9 


15.2 


2.7 


208pb 


RPA 


12.7 


12.9 


13.4 


1.8 
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TABLE II: Mean energies and widths for the ISGMR strength distributions in the even-A Sn 
isotopes calculated within the QRPA for 10.5-20.5 MeV energy interval. The calculations were 
performed within self-consistent HF+BCS approach based on the T6 Skyrme force (/Coo = 236 
MeV). The RPA results for ^°°'^^^Sn and 208p^ g^^^ g^jg^ gj-^own. The values for 208p^ g^^^ calculated 
for 5~25 MeV energy interval. 





-\ / 771 ^ / 771 i A/TpV 


771 1 1 711 r\ A/TpV 


•\ / 777 o 1 771 1 A/TpV 


r, MeV 


lOOc 

kJli 


1 7 4 


1 7 5 


1 7 8 


9 4 


"2Sn 


17.1 


17.2 


17.4 


2.2 


ii^Sn 


17.0 


17.1 


17.3 


2.3 


"6Sn 


16.8 


16.9 


17.2 


2.4 


"8Sn 


16.7 


16.8 


17.1 


2.5 


i20Sn 


16.6 


16.7 


17.0 


2.7 


i22Sn 


16.5 


16.6 


16.9 


2.8 


i24Sn 


16.3 


16.5 


16.8 


3.0 




16.0 


16.1 


16.5 


3.2 


208p|r, 


13.9 


14.1 


14.7 


2.0 



As can be seen from Table UTTl the (A^ — Z) dependence of the (Q)RPA mean energies prac- 
tically follows the dependence of the IQM energies. In particular, the difference between the 
mi/rriQ values for ^^^Sn and ^^^Sn in the QRPA is equal to 0.9 MeV and the same difference 
is obtained in the IQM calculation. Since the poles of the uncorrelated propagator A{uj) 
are equal to the sums of the quasiparticle energies + E(^2) [see Eq. (IA3P ]. this result 
means that the (A^ — Z) dependence of the ISGMR mean energies is mainly determined by 
the level density of the single-quasiparticle spectrum. Including the residual interaction in 
the (Q)RPA, we obtain the following redistribution of the isoscalar monopole strength: the 
low-lying part of the strength disappears, while the mean energy of the high-lying states 
(which form the ISGMR) is reduced by approximately two MeV. 

In Table UTTj we also include the ISGMR mean energies obtained within the RPA for ^^°Sn 
nucleus. In this calculation, the pairing correlations are neglected both in the mean field and 
in the residual interaction. Respective strength function is shown in Fig. [2] in comparison 
with the IQM and QRPA strength functions. These results demonstrate that the influence 
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TABLE III: Mean energies for the 0^ strength distributions in the even- 

^ 100,112-124,132 isotopes 

calculated for 10-30 MeV energy interval within self-consistent HF-I-BCS approach based on the 
T5 Skyrme force. See text for details. 







a/?7?,i /m, 1 MeV 


777,1 1 Tn.r\ A/TpV 


\ / Tn,i 1 ?7?,i MeV 


lOOSn 


IQM 


19.1 


19.3 


19.9 




RPA 


16.8 


17.0 


17.6 


ii^Sn 


IQM 


18.5 


18.7 


19.4 




ORPA 


16.4 


16.5 


17.1 


ii^Sn 


IQM 


18.3 


18.5 


19.3 




QRPA 


16.2 


16.4 


17.0 


i^^Sn 


IQM 


18.1 


18.4 


19.1 




QRPA 


16.0 


16.2 


16.8 


ii^Sn 


IQM 


17.9 


18.2 


19.0 




QRPA 


15.9 


16.0 


16.7 




IQM 


17.8 


18.1 


18.9 




QRPA 


15.7 


15.9 


16.5 




RPA 


15.2 


15.4 


IC.O 


i22Sn 


IQM 


17.6 


17.9 


18.7 




QRPA 


15.6 


15.7 


16.4 


i24Sn 


IQM 


17.5 


17.8 


18.6 




QRPA 


15.4 


15.6 


16.3 


i32Sn 


IQM 


17.1 


17.4 


18.2 




RPA 


15.0 


15.2 


15.8 



of the pairing correlations on the ISGMR mean energies is appreciable. In the QRPA, where 
the pairing correlations are included, the mean energies increase by 0.5 MeV as compared 
with the RPA. By comparing the results obtained with the T5 and T6 Skyrme forces, one 
can see that this increase is substantial for determining the INM incompressibility. 
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Energy (MeV) 

FIG. 2: Isoscalar EO response in ^■^''Sn calculated within the independent quasiparticle model (IQM, 
dashed line), RPA (dashed-dotted line), and QRPA (solid line), making use of the T5 Skyrme force. 
See text for details. The smearing parameter A is equal to 500 keV. 

IV. CONCLUSIONS 

In the paper the results of the theoretical analysis of the ISGMR strength distributions 
in the even-A ^^^~^^^Sn isotopes are presented. The calculations were performed within 
two microscopic models: the quasiparticle random phase approximation (QRPA) and the 
quasiparticle time blocking approximation (QTBA) which is an extension of the QRPA 
including quasiparticle-phonon coupling. We used self-consistent calculational scheme based 
on the HF+BCS approximation. The self-consistent mean field and the effective interaction 
were derived from the Skyrme energy functional. In the calculations, two Skyrme force 
parametrizations were used. The T5 parametrization with comparatively low value of the 
incompressibility of infinite nuclear matter {Koo = 202 MeV) allowed us to achieve good 
agreement with the experimental data for tin isotopes within the QTBA including resonance 
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widths. However, this parametrization fails to reproduce the experimental ISGMR energy 
for the ^°^Pb nucleus which is usually used in the fit of the Skyrme force parameters. On 
the other hand, the T6 Skyrme force with = 236 MeV nicely reproduces the ISGMR 
energy for 208p|^ overestimates the energies for ^^^~^^^Sn isotopes by more than one 
MeV. On the whole, these results do not allow us to decrease the ambiguity in the value 
of Koo as compared with the previous known estimates. Note, however, that the main goal 
of our work is not to solve the problem of the nuclear matter incompressibility but to find 
under which conditions one can obtain reasonable description of the experimental data for 
the considered tin isotopes within the framework of the self-consistent approach including 
correlations beyond the QRPA. 
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APPENDIX A: MODIFICATION OF THE QTBA EQUATIONS INCLUDING 
CONTRIBUTION OF THE PARTICLE-PARTICLE CHANNEL IN TERMS OF 
THE REDUCED MATRIX ELEMENTS 



In the detailed form using the notations of Ref. [23|] for the reduced matrix elements, our 
method to include pp-channel contribution in the QTBA equations consists in the following. 
In Eq. (33) of Ref. [231 only ph channel is kept, but in Eq. (42) for Aj^^^^J^P^^^^'^'^V) the 
matrix element ylj^2 34]('^) is replaced by A^^''^^r^^'^^\uj) where 



.J(res+pp)/ X 

^[12,34] y^) 



A A /(■^(res+pp) / \ 

% .-'72 %-V4 ^(12)r,,, (34)773 ) 



(Al) 



Propagator ^(i27,^^34)V ^ solution of the equation 



A J (res+pp) 
^(12)^, (34)7, 



^(12)7,, (34)7?' 



-J (res+pp) 



{u;)A 



J (res+pp) 



+ ^ ^^(65) ^(12)7j, (56)7?" ^■'(56)7?",(34)7j'' 



(A2) 



7?" (56) 
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where 



(fsJ (res+pp) / N 



J (res+pp) 



■'?l'?2'?3'?4 



(A5) 



J (res+pp) 



[12,34] V^J '*'[12,34] V"^7 [12,34] [12,34] ' 

The order-bounding factors O^^^i) Eq. ( ]A2l) are defined as follows: 0^.2]^) 
number of the state (1) is lesser than the number of (2) [(1) < (2)], 9,^-,^ = ^ if fl 

e 



^[12,34] l^J 



(res) /p,N W(PP) 
^[12,34] l*^^ "t" -^[12,: 



'(21) 



(A6) 

1 if the ordinal 

(2), 



'(21) 



if (1) > (2). The interaction amplitude ^'^^^m]^^) responsible for the QPC is 



defined by Eq. (B14) of Ref. [23|. Introducing the notation 



(12)r,, (34)r,' 



V1V2V3V4, 



J(pp) 



»7j,r) r?2,-r? i)^,r]' r?^ [12,34] 



(A7) 



and using Eqs. (C2)-(C4) of Ref. [23| we obtain the following ansatz for this quantity 



r, 



j(pp) 



-02^2||Tjjo ||jl/l)(j4^4|| Tjjo ||j3^3) 



(12)r?, (34)r?' , "^g , ^4 , ^3 2 J + 1 

^ \^r],'r]' (^(1)^(2)'"(3)^(4) + "^(1)^(2) "^(3) ^(4)) ~ ^'r),-r)' (""(1)^(2) "^(3) ^(4) + '^(1)'^(2) ^(3)'"(4) )] 
/■oo 

X / c?rr2i?(^^(r)/?(2)W^(3)W^(4)W-^^W- (A8) 

t/ 

Note that the value of -^(12)^^ (34),,' Eq. (lASP of the present paper differs from the 



corresponding value derived from Eqs. (C2) and (C3) of Ref. |23| by a factor \ due to the 
shorthand summation used in Eq. (CI) [(3) ^ (4)]. In addition, in the case J = one 
should set: ^(12)!^^ (34)^/ =- '^(12) '^(34) -^{xiy^ {zz)ri' order to obtain consistency with the gap 
equation (A25) of Ref. [23i] written in the diagonal approximation. Note that this method 
is applicable both in the QTBA and in the QRPA. In the latter case the amplitudes (^^^"^^^^ 
in Eq. flA6P are set to be equal to zero. 
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